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Abstract 

We investigate the description of statistical field theories using Daubechies' 
orthonormal compact wavelets on a lattice. A simple variational approach 
is used to extend mean field theory and make predictions for the fluctuation 
strengths of wavelet coefficients and thus for the correlation function. The re- 
sults are compared to Monte Carlo simulations. We find that wavelets provide 
a reasonable description of critical phenomena with only a small number of 
variational parameters. This lets us hope for an implementation of the renor- 
malization group in wavelet space. 



1 Introduction 



The success of wavelets [1-6] in analyzing complex signals has prompted speculation 
about their possible applications to field theories and other lattice systems. Like 
the Fourier transform, the main feature of the wavelet transform is the separation 
of scales. It keeps, however, a (coarse-grained) position-space lattice that enables 
one to mix the features of real- and Fourier-space block-spin transformations, and 
may thus be a possible candidate for a renormalization group that could be more 
efficient both in terms of how accurate physical processes are depicted and how 
much numerical work is necessary to calculate the renormalization group flow. The 
localization feature of wavelets is especially appealing when one considers theories 
that have a long-range scale (as the confinement scale in Quantum Chromodynamics). 

*email: cbestOth. physik. uni-f rankfurt .de 
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To find out whether a renormalization group description would work in wavelet space, 
one must learn about how efficiently a field theory can be described in terms of 
wavelets. 

In this paper, we study the representation of statistical field theories in wavelet 
space on the level of the Gaussian approximation. We use a simple variational ansatz 
to calculate fluctuations in wavelet space. These are related to the correlation func- 
tion and can be compared to results of Monte Carlo calculations. Other variational 
approaches based on the Hamiltonian formalism |7j have been quite successful in 
Quantum Chromodynamics (QCD). 

In section II, we give a short introduction to the wavelet transform. Our varia- 
tional procedure is defined in section III. Its application to the free field theory and 
to the XY model is demonstrated in section IV. In section V, we apply the method 
to a Landau- Ginzburg field theory. 



2 The wavelet transform 
2.1 Definition 

A wavelet ip G L 2 {R) is a function whose binary dilatations and dyadic translations 
generate a Riesz basis of L 2 (R). That is, any function in £ 2 (R) can be expanded 
into a wavelet series, 

oo 

/(*)= E E / (n V)v> (n V)(*) , (i) 

n=— oo x'£C" 

where ip^ n \x') £ L 2 (R) denotes the dilatated and translated wavelet, defined by 

ijjW(x')(x) = 2~ n ' 2 ij (2- n {x - x'j) . (2) 

iiGZ gives the scale of the wavelet, corresponding to a dilatation by 2 n , and x' E C n 
the position (translation) on the sublattice C n = 2 n Z of scale n. 

The coefficients p n \x') characterize features of the wavelet at scale n and po- 
sition x' . In this sense, the wavelet transform offers a mixed position-frequency 
representation of the signal. 

Compact orthonormal wavelets have been constructed by Daubechies M. They 
cannot be given in closed form but by a powerful numerical procedure || [Hj. The 
wavelet transform constitutes a multiresolution analysis in which the function space 
L 2 (R) is decomposed into a set of nested subspaces 

• • • C Vi C V C V-x c • • • (3) 

generated by binary dilatations of a scaling function <p(x): 

V n = dos [<p {n \x') : x' e 2 n Z} (4) 
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Besides the requirement that their union covers the whole function space and their 
intersection is the zero function, they have the important property that the contrac- 
tion f(2x) of any function / G V n is a member of the next finer subspace V n+ \. As 
it turns out, the pair of scaling function and wavelet can be chosen such that the 
spaces V n are direct sums of (disjoint) wavelet spaces, 

W n = clos{ilJ {n) (x , ):x , e2 n z} , (5) 

V n = V n+1 ®W n . (6) 

Wavelets thus generate a decomposition of the function space into disjoint subspaces 
W n , representing different scales. Unlike other methods of multi-resolution analysis, 
the wavelet spaces and their scaling function spaces are derived from a single pair of 
functions, ip{x) and (j)(x), by simple translations and dilatations. 

2.2 Construction 

The construction of wavelets proceeds from the property that the scaling function 
(f)( n ) G V n and the wavelet ift( n ' G W n are both members of the next finer subspace 
V n -\- Thus, they can be expressed as a combination of the basis function {0^ n " 1 ^(x / ) : 
x' G C n - 1 } of V n - X : 

^ n \k){x) = ^/i i (n - 1) (A; + 2"- 1 /)(a;) 
i 

^ n \k)(x) = Y.di^ik + T-H^x) (7) 

This is the two-scale relation. The coefficient series {/iz}z 6 z and {gi}iez characterize 
the wavelet. Compactly supported wavelets have a finite number L of nonvanishing 
coefficients. For orthogonality, the coefficient series must be related by 

g l = (-\) l h L -i-x ■ (8) 
Eq. (|7j) can be inverted to give [x' G £ n_1 = 2 n_1 Z) 

^ n ~ 1 \2x')(x) = Y.iS^W ' -2 n l){x)+g 2l ^ n \2x' -2 n l)(x)) 



i 



(n-l)/o„/ 



(2x' + l)(x) = J £{h a+1 ^ n \2x'-2 n l){x) 



i 



+g2i+i^ n \2x' -2 n l){x)) . (9) 



Eq. (0) is the basis of the partial wavelet transform. If we consider the scalar 
product of a function f(x) with scaling functions and wavelets, resp., 



f {n \x') = J (f) {n \x){x)f{x)dx , (10) 
f {n \x') = f^ in) {x'){x)f{x)dx , (11) 
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the coefficients f( n \x') give a smoothed representation of f(x), while the wavelet 
coefficients p n \x') contain the complementary detail information that is lost in the 
smoothing process. Using eq. (0), both can be expressed by smoothed coefficients on 
the next finer scale, 

/ (n V) = 52h l F n - 1 \J + 2 n -H) , (12) 
i 

/ (n V) = Y,9if {n - 1 \x' + 2 n - 1 l) . (13) 
i 

Repeating the procedure, any wavelet coefficient can be reduced to the smoothed 
coefficients at as fine a scale as desired, essentially corresponding to the resolution of 
the measurement. Using the inverse transformation @, the basis functions (wavelets 
and scaling functions) can be reconstructed as accurately as desired. 
Eq. ([12]) defines an orthogonal transformation between the set 

{/ (n_1) 0O = J e ^} (14) 

of smoothed coefficients at scale n — 1 and the set of smoothed plus detail coefficients 
at scale n 

{f {n) (x')J^(x'):x' EC n - 1 } . (15) 

If we consider a finite-dimensional vector (f(°>(x'),x' — 0, . . . , 2 N — 1) as input (imply- 
ing periodic continuation where required), repeated application of this transformation 
results in the wavelet vector 

{f^(x'),n = 1, . . . , N, x' G C n ; f {N) (0)} (16) 

As we have only a finite number of elements in the input vector, the partial wavelet 
transform must stop when there is only one coefficient left, the smoothed coefficient 
of the coarsest scale, giving the average of f(x). The transformation is orthogonal 
and thus defines an orthonormal basis, consisting of discretized wavelets, in the space 
of vectors. The scaling relation (0) applies only approximately to these wavelets. The 
higher the scale n and the larger the extent of the wavelet, the less is the discretization 
felt, and the discrete wavelets approach the continuum ones. 

Wavelets in more than one dimension can be constructed by direct products. They 
then carry an additional label t — 1, . . . , 2 D — 1 that specifies their composition. It 
is advantageous to put both wavelet and scaling function on the same footing using 
an index t — 0,1, 



^ j ^W( X ')( X ) , if t = 1 (17) 



Then the iV-dimensional scaling function (t = 0) and wavelets (t = 1, 
are given by 

4 n \x')(x) = f[4t D ' n \x' k )(x k ), 

k=l 



•>D 
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where 



t = E t k 2 k -i 

k=i 



(19) 



3 Method 

3.1 Wavelet representation 

We study statistical field theories on the lattice in wavelet representations. The 
lattice field a(x), x G C, is expanded over a wavelet basis, 

N n w 

«(*) = EE E at\x')4 n \x'){x) (20) 

n=l « x'e£ n 

^V) = E^MW»W > ( 21 ) 

on a D-dimensional lattice £ = £° with 2 DAr sites, A?" giving the number of scales. 
The sublattices C n have accordingly 2 D ( N ~ n ^ sites, t = or 1, ...,n w = 2 D — 1 
denotes the cartesian composition of multidimensional wavelets, t = is included 
only at the topmost level n — N, where (f)^ is the constant function, and thus 
gives the average of a(x) over all sites. 

3.2 Variational principle 

We employ the principle of minimal free energy to obtain an approximate description 
of the probability distribution in the partition sum. In a trial ensemble with prob- 
ability distribution P(a; {a(x)}) characterized by a set of variational parameters a, 
we calculate the entropy S 

S = J d{a(x)}P(a,{a(x)})\nP(a,{a(x)}) (22) 

and the internal energy U 

U = J d{a{x)}P(a,{a{x)})H{{a(x)}) . (23) 

By minimizing the free energy 

F = U-±S (24) 

with respect to the parameters a, the best-fit probability distribution is obtained. 

As our goal is to find out how efficiently wavelets describe statistical field theo- 
ries, we choose as a trial ensemble the coarsest possibility, a Gaussian uncorrelated 
ensemble in wavelet space, 

(at\ Xl ) at\x 2 )) = 5 nun2 5 tut2 8 XUX2 (25) 
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with scale-dependent fluctuation strengths A\ ■ This ensemble makes calculations 
sufficiently simple and is, at the same time, able to describe nonlocal correlations in 
position space, as we show below. It can be thought as an extension of the purely 
local ansatz 

(a(x)a(y)) = 6 x , y A (26) 

used in the simplest Gaussian approximation, which corresponds in wavelet space 
to constraining all A^ to the same value. The introduction of different fluctuation 
scales in the wavelet transform hence enables us to improve the description of nonlocal 
fluctuations. 

The entropy of this ensemble is easily calculated to be 

i N n w 

S-^^iVjEi!" 1 + const . (27) 

^ n=l t 

N n = \C n \ = 2 D( - N °~ n > is the number of sites on the n-th sublattice. 
3.3 Correlator 

Given the fluctuation strengths A[ n \ the correlator can be calculated by applying 
the wavelet transform 

(«W«(!/))=EE^ ) (^)A (n) (28) 

n t 

with the elementary wavelet correlator 

Ct n) (x,y)= £ 4 n \x')(x)4 n \x')(y) • (29) 

These functions are not translation invariant as the translation group has a com- 
plicated representation in wavelet space. Translation invariance of the correlator 
(f28|) must thus be restored approximately by the A^ in the same way as rotational 
invariance is restored on a lattice correlation function. 

The scaling relations between wavelets (which hold approximately for discrete 
wavelets of sufficiently high scale n) imply 

Ci n \x, y) « d n \z - y) ~ 2^«)/ 2 (2 n '~ n (x - yj) . (30) 

The elementary wavelet correlators have finite but nonzero extent (as the wavelets 
have). Thus, long-range correlations can be obtained by (in wavelet space uncorre- 
lated) fluctuations at sufficiently large scales. In this way, wavelets can map a critical 
position-space theory with long-range correlations to a less critical wavelet-space the- 
ory dominated by short-range correlations. 
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4 Gaussian and XY model 



4.1 Wavelet fluctuation strengths 

The Gaussian model (free field theory) is given by the Hamiltonian 

w = ^EE(^ + ^-«W) + y^( a(x)2 ) • (31) 

x&C n=l x£C 



It can be solved exactly by a Fourier transform, giving the correlator [ ] I 



1 1 x—* cos(/c ■ d) 

«»W* + <*)> = E £ ml + 2 j^ (2 _ 2 cos(fc . ^ (32) 

with fc on the dual lattice JF£ with lattice spacing 27r/2 N . 

The wavelet correlator can be calculated directly by applying the wavelet trans- 
form to this expression: 

(a^Vi)4" 2) (4)>= E ^ ( r ) K)(x 1 )^ 2) (4)(x 2 )(a(a; 1 )a(x 2 )) . (33) 

X1,X2&C 

On the other hand, fluctuations in wavelet coefficients can be obtained from the 
variational principle. To this end, we use ( P5| ) and calculate the internal energy in 
the ensemble. To simplify the calculation, we rewrite the next-neighbor difference as 
a scalar product in lattice field space between the field and the characteristic function 
of the link, 

a(x + e M ) - a(x) = (a, L(x, /i)) = a (v) L ( x ' ^(v) ( 34 ) 
with the characteristic function 

L (X,V)(V) = <We„ - 6 Via . (35) 

As the wavelet transform is orthogonal, it conserves the scalar product 

N n w 

a(x + — a(x) = (a, L(x, //)) = E E E a i n \ x ') Lt\x, ^){x') (36) 

n=l t x'eC n 

where L[ n \x, n)(x') is the wavelet transform of L(x,fi)(y). Then the expectation 
value of the link Hamiltonian is 

i \ i N n w 

-(a(x + e,)-a(x)f) = - £ £ $»W)M b) (37) 
with the symbolic notation 

|L^(x,/i)| 2 := £ $ B W)(*0) a (38) 
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This gives the power spectrum of the link (x, fi) on scale n, in the sense that the 
original norm of the function L(x,fi) is distributed over different scales, subject to 
conservation of the norm 

N n w 

££|Li B W)l 2 = E(£(*>A0(i/)) a ■ (39) 

n=l t yeC 

The notation can be extended to 

:=££|^W)| 2 (40) 

which, so to speak, gives the power spectrum of the whole lattice. 
Similarly, the site Hamiltonian gives 

m 2 m 2 N nm 

^) 2 = ^EEi^ (n) (-)r^ n) (4i) 

Z Z n=l t 

with the characteristic function S(x) of a site: 

S(x)(y) = 6 VtX . (42) 

In this case, there is further analytical simplification as the sum over x G C can be 
performed: 

E (Sl n \x)(x')) 2 = E (#V)(*)) 2 = 1 (43) 

due to normalization of the wavelets. 
The internal energy is therefore 

1 N 

^ = Wo = ^EE^ n) l 2 ^ n) + ! fEE^^ n) • (44) 

/ n=l t z n=l t 

Both entropy and internal energy are linear in the variational parameters, so that 
the variation can be performed analytically, leading to the result 

4 n) = \ ^1 (45) 

with the notation ([TOp. The quantity \L[ n ^\ 2 /N n has taken the place of the momentum 
square k 2 in the corresponding Fourier-space expression. 

The wavelet coefficient n = N, t = gives the average of the function over all 
sites. In the case of a characteristic function of a link ( |35|) , the average vanishes and 
thus 

4 W) = ^ (46) 
8 



Gaussian correlator, m = 0.032, 64 x 64 

~i 1 1 1 1 r 

exact — 

DAUB4 --- 

DAUB8 

DAUB20 --- 




5 10 15 20 25 30 

d 

Figure 1: Correlator of a free field at f3 — 1. The solid line is the exact result, the 
broken lines the result of the variational procedure with different wavelet types. 

is divergent for m? — > 0. The same divergence appears in the Fourier transform 
associated with the k = 0-mode and is caused by the invariance under adding a 
constant to a massless field. As it appears only in a single wavelet coefficient, it can 
be as easily subtracted in the wavelet with the Fourier transform. 

Fig. m shows the prediction for the correlator from different wavelets. The higher 
the order of the wavelet, the more correct the result. This is not surprising as 
higher wavelets approach more and more the form of free waves that diagonalize the 
Hamiltonian. 



4.2 Correlator 

To calculate the correlator of the massless theory, we first derive the scaling properties 
of the fluctuation strengths At in respect to n. Using the factorization property 
(0) of the wavelets, we find for (gDJ) 

Lt\x^){x') = 4 n \ x + eix -x')-4 n \x-x') 

x(4i D ' n \x» + l-x'»)-4l D ' n) M) (47) 

Summing over x', x, and /i yields 

l^i 2 = 2 D ^Dj2(4t D ' n \x, + 1) - 4l D ' n) M) 2 . (48) 
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The one-dimensional finite difference can be approximated by a derivative as the 
discrete wavelets approach the continuum wavelets (@) for large n 

4i D ' n \x» + i)-4t D ' n \x») 

« 2- n/2 (^(2" n (x + 1)) - ^j(2- n x)) » 2- 3n/ V(2" n x) . (49) 



Consequently 



X 

« 2" 3n J] (^(2~ n x)) 2 « 2~ 2n f dx (^'(x)) 2 (50) 

a: ^ 

and we deduce the scaling relation 

An 

|£(n)|2_ 2 -(D+2)n ,(n) ^ 4 

i t i * 2/3 v ; 

We also note here that, in the massive case m^O, ( fl5|) implies that the .4^ tend 
to the constant value l//3m. As we shall see, this corresponds to the absence of 
long-range correlations. 

For the wavelet autocorrelation c[ n \x,y) function, we capture its gross features 

by 

C^(x,x + d) « (1 - 2d) 9(1 - d) 9(d) . (52) 
The correlation function of the field theory is then 

AT n-u, 

(a(x)a(x + d)) = $]Z!^S n) (l-2 1 " n d)e(l-2 1 - n d)e(d) . (53) 



n=l t 



(n) 



Let 2 n ° < d < 2" 0+1 . Then only terms n > uq contribute in the sum. If A, 
a n A^°\ the sum can be performed, sending iV — > oo, 

oo 

C(d) = A {0) n w e(d) a n 2- Dn (l -2 l ~ n d) 

n=no 

= A^n w Q(d)a^2-^ _ 2 i-«o d _^_^ . (54 ) 

Putting 2 no ~ d and thus a n ° ~ d log2 ", we find 

C(d) ~ d" (1) - 2) . (55) 

This is the expected scaling behavior for a Gaussian theory [II | . In two dimensions, 
it corresponds to a logarithmic correlation function. 
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4.3 XY model 



We would like to test our approximation on a nontrivial field theory. The two- 
dimensional XY model is among the simplest field theories with a nontrivial critical 
behavior. 

The Hamiltonian of the XY model is 

D 

Tt = cos (a(x + e M ) — a(x)) . (56) 

The field — tt < a(x) < tt gives the angle of the spin at x, relative to some global 
reference direction. It is symmetric under global 0(1) transformations and undergoes 
a Kosterlitz-Thouless phase transition in two dimensions [ |ll|] . 

As in the Gaussian model, we expand the field a(x) in a wavelet basis, assume 
Gaussian local fluctuations, and calculate the internal energy as the expectation 
value of the Hamiltonian in the ensemble. However, to make a Gaussian ansatz, 
we have to extend the field a(x) to the whole real line. The Gaussian ansatz now 
does not reflect the periodicity of the Hamiltonian. This is a valid approximation in 
the ordered low-temperature phase, but it fails at high temperatures. This fact is, 
however, independent of the wavelet expansion. 

To calculate the internal energy, it is useful to note that 

(cos(x + x )) = exp I — J cosx (57) 

for a Gaussian variable x. With this and the same procedure as above, we find for 
the internal energy 



^ = EE 



/ i N n w n 



(5? 



In the limit of small fluctuations, the exponential can be expanded, and the (mass- 
less) Gaussian model is recovered for the low-temperature limit. Conversely, at high 
temperatures the exponential vanishes and the internal energy approaches unity per 
site, reflecting maximal disorder. 

The free energy exhibits in this approximation a local minimum that vanishes 
when the temperature increases. The absence of a global minimum is rooted in the 
lack of periodicity of the Gaussian probability distribution. While the internal energy 



U is bounded by unity at maximal disorder, our entropy (|27j ) can increase without 
bounds. The true entropy, however, is bounded by the available configuration space, 
'S'max — hi 27r. By assuming a Gaussian probability distribution, we have extended 
the configuration space to the infinite real line. 

We can, however, still interpret this result from the point of view of ergodicity 
in a system with a periodic potential: In the low-temperature phase, the system is 
localized in a small part of the available configuration space. When the temperature 
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Figure 2: Internal energy per site of the two-dimensional XY model on 64 x 64 
lattice, calculated variationally with wavelets and by a Monte Carlo method. 



increases, the localized solution becomes unstable at some point, but as our trial 
probability distribution space does not include periodic distributions, the minimiza- 
tion problem has no solution. 

Numerically, we find the temperature of this phase transition at /3 cr it ~ 0.69 in 
two dimensions and (3 CT \ t ~ 0.45 in three dimensions. This is consistent both with 
the result for the purely local ansatz (|26|) , which gives /3 crit = e/2D « 0.6769 ... in 
two dimensions and 0.4530 ... in three dimensions, and with the prediction from an 
analysis of vortices in two dimensions giving /? crit = 2/n « 0.6366 . . .. 

Wavelets, however, do not only predict the internal energy but also nonlocal 
correlations that are not included in the position-space ansatz. Fig. [$] compares the 
lowest wavelet fluctuation strengths as calculated variationally and as measured in 
a Monte Carlo simulation. This essentially amounts to comparing correlators, which 
are not predicted by In the Monte Carlo simulation, we have ignored periodicity 
(i.e. allowed a(x) to be outside the interval [— 7r, tt]) to model the Gaussian ansatz as 
closely as possible. 

Above (3 m 1.5, variational and Monte Carlo result coincide. Below that, nongaus- 
sian fluctuations due to periodicity begin to play a role, and the variational method 
fails. This is also manifest when ([a[ n \x')] A ) is measured where a strong departure 
from Gaussian factorization can be observed. 
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Wavelet fluctuation strengths, XY model, 64 x 64 (angles) 
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Figure 3: Wavelet coefficient fluctuations in the two-dimensional XY model on a 
64 x 64 lattice. 



5 Landau-Ginzburg model 

In the example of the XY model, major problems were caused by having the (peri- 
odic) group parameter (i.e., the angle) as field variable instead of the spin vectors. 
Since the constraint s 2 = 1 on the length of the spin vectors translates to a com- 
plicated expression in wavelet space, it is not easily incorporated in the variational 
ansatz. (The same problem prevented us from considering the even more simple Ising 
model.) 

However, as our main interest will be in renormalization group transformations 
which dilute the constraints and eventually allow a description in terms of a contin- 
uous order parameter governed by an effective theory, we now turn to the study of 
this effective theory, the Landau-Ginzburg model. 



5.1 Wavelet representation 

The Landau-Ginzburg model is given by a iV c -component field a(x) with the 0(N C )- 
symmetric Hamiltonian 

in 2 

w = .EE («(* + - «(*)) + \ E + ?E («(*) 2 ) (59) 

z xec^=i z xec z xec 

corresponding to a ($ 2 ) 2 field theory. For m 2 < and g > it exhibits a "Mexi- 
can hat" potential with spontaneous symmetry breaking. In mean-field theory, the 
phase transition between 0(N C ) symmetric and symmetry-broken phase takes place 
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at (3 — 0. We wish to estimate the shift of the phase transition temperature due to 
fluctuations in wavelets. 

To allow for symmetry breaking, we must add another variational parameter, the 
expectation value a of the field: 

(a(x)) = a . (60) 
As this is a constant, it implies for the wavelet coefficients 

(at\x'))=5 n , N 8 t , Q -?= , (61) 
V iV o 

i.e., the only coefficient with a nonzero expectation value is clq N \ associated with the 
scaling function at the topmost level. The additional variational parameter a, giving 
the center of the probability distribution of ajj , appears only in the internal energy, 
as a mere shift of the probability distribution does not influence the entropy. 

The first two terms of the Hamiltonian (kinetic and mass) have been calculated 
in the context of the Gaussian model. The only modification is that we now have to 
deal with a vector- valued field a a (x), a = 1, . . . , N c . The wavelet fluctuation strength 
is thus a tensor 

(^l) ^is^aO*^)) = ^ni,n 2 ^ii,i2^a;' 1 ,x^ • / ^i 1 ^a 1 a 2 • (62) 

The unperturbed Hamiltonian appears then as 

1 N n w 2 N n w 2 

Ho= EE l^ n) I 2 tr ^ + ^-EE Ar « tr 4° + a 2 (63) 

1 n=l t 1 n=l t 1 

where the trace operator refers to summing over the vector indices a. 

To calculate the influence of interactions, we expand the third term of ( |59| ) in 
wavelet space and then can make use of Wick's theorem, as our trial ensemble is 
Gaussian: 



Then 



(abed) = (ab) (cd) + (ac) (bd) + (ad) (be) . (64) 
E ( E a a( x ) 2 ) = E E aa(x)a a (x)ap(x)ap(x) 

x \ a ) x a/3 

z ni,n 2 tx,*2 x 

x(tr^ 1 1 tr^ + 2tr(ATA n 2 2 )) 

N n w „ 

+s E E E E K n) (*') (*)] « 2 A (n) + 2 a ( ;A) 



n t x' % 



+ 9 -N a A . (65) 
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We take the symmetry-breaking along the 1-direction in 0(N C ) space. The fluctuation 
matrix at each scale retains the symmetry under the corresponding stabilizer group 
(i.e. the group of transformations that leave the symmetry breaking vector invariant) 
and is thus of the form 



A 



(A ••• \ 

A± ••• 

o o Al ••• 

V ; ; i ■■ J 



(66) 



A\\ gives the fluctuation strength in the direction of the symmetry breaking, Aj_ 
perpendicular to it. In particular 



trA = A\\ + (N C -1)A± 
tr(AA') = A\\A\ + (N c - 1)A±A' ± 

With this, the complete internal energy is 

1 JV n„, 

U = ^E(\ L ^\ 2 + m2N n)(^\S + (N c -l)A 



(67) 
(68) 



(«) 



n=l t 

N n u 



+ | E E MMS + ^ ~ ^SMS 



|(«i) /i(n 2 ) 



n\U2 ti*2 



N n u 



+2(N C - 1)A^A^ 



(ni) A (n 2 ) 



n t 



^ a 4 



(69) 



The first term gives the original Gaussian Hamiltonian, the second describes inter- 
actions between different scales due to the $ 4 interaction term; the final two terms 
are due to spontaneous symmetry breaking. 
We have introduced the matrix 



M nin 2 

t\t2 



E E E KrViX 



x: 



2 r 



(70) 



describing the nonlinear interaction of fluctuations at different scales caused by the 
quartic term. Due to the factorization properties of wavelets, it can be rewritten as 



D 



n m% 



lfc*2fc 



k=l 



with 



(71) 



(72) 



x' 2 £C n 2 x£C 



derived from one-dimensional wavelets. In particular, the complexity of (|69|) increases 
only slightly when higher dimensions are considered. 
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5.2 Spontaneous symmetry breaking 



The value of a can be calculated analytically. As it does not enter the entropy, it 
is determined by the minimum of the internal energy if the A[ n) are given. The 
corresponding condition is 



dF 2 N n 



in-Noa + 2N ga 3 + 2ga £ £ N n (3A$ + (N c - l)A^ t ) = . (73) 



da „ , 



While a = is always one solution to this equation, the other one is the root of a 
quadratic equation, 



m 2 N nw , , \ i \ 

^-^-EE^(3^ + (iv c -i)<; 



(74) 



Without the fluctuations, this is just the mean field result which exhibits sponta- 
neous symmetry breaking for all values of (3, provided that m 2 /2g < 0. The more 
fluctuations appear, the more a is reduced, until the expression under the root be- 
comes negative, and a = becomes the only solution, corresponding to a symmetric 
phase. 

The value of a 2 interacts with fluctuations due to the third term in (|69|). As 
the fluctuations enter linearly and indiscriminately of their scale, this modifies the 
effective mass of the theory. Due to the prefactor 3 in the fluctuation sum, it is 
different for parallel and for perpendicular fluctuations. Inserting the mean-field 
value of a for comparison, we find 

- m + Qga 2 w —2m , , 

= m + 2ga 2 w (75j 




The value —2m is the correct Gaussian curvature of the potential for fluctuations 
against the valley. Fluctuations perpendicular to the magnetization are massless 
Goldstone modes. 



5.3 Numerical solution 

We tested the accuracy of this model both for the case of a small anharmonic per- 
turbation and the full phase transition. Eq. (]69D was minimized numerically using a 
combination of the simplex method and simulated annealing ||. 

For comparison, a simple Metropolis simulation of the model was performed. The 
wavelet transform was calculated on statistically independent configurations and the 
fluctuation strength of wavelet coefficients measured in this ensemble. 

One should keep in mind that the trial probability distribution used here makes 
at least two bold approximations: It neglects all correlations between wavelet fluctu- 
ations, and it assumes that all fluctuations are Gaussian. We therefore should expect 
no more than a gross representation of the model. However, this should be sufficient 
for a wavelet renormalization group to work, as the renormalization transformations 
should enhance the physically important features. 
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Wavelet coefficients, anharmonic perturbation 




5.3.1 Small anharmonic perturbations 

Small values of g > represent anharmonicities that manifest themselves in the 
departures from the free-field result, comparable to perturbation theory. Fig. f| shows 
the wavelet coefficient fluctuations calculated variationally and compared to the result 
of a Monte Carlo calculation. Even where the departure from the free-field result is 
large (at small (3), the results are close. 

5.3.2 Landau-Ginzburg phase transition 

For g > 0, m < 0, the potential has the "Mexican hat" potential leading to sponta- 
neous symmetry breaking in mean field theory. Fluctuations change the temperature 
of the symmetry-breaking phase transition from the mean-field value /3 — to some 
finite value. In Fig. [5], we have calculated the magnetization variationally both using 
the wavelet basis and the localized ansatz (|26| ) and by a Monte Carlo calculation on 
a 32 x 32 lattice with parameters m? = —4 and g = 2. In the Monte Carlo calcula- 
tion, an external field h = 0.05 was applied. The finite lattice size makes the phase 
transition less pronounced, but the inclusion of nonlocal fluctuations in the wavelet 
ansatz lowers, as expected, the phase transition temperature towards the result of 
the Monte Carlo calculation. 

Fig. P shows the wavelet fluctuations calculated variationally and in the Monte 
Carlo simulation. The variational wavelet ansatz gives a faithful representation of 
the overall features. In particular, the wavelet coefficients both in the Monte Carlo 
and the variational calculation develop an approximate power scaling near the phase 
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Figure 5: Magnetization in the Landau- Ginzburg model on a 32 x 32 lattice with 
m 2 = —4, g = 2, calculated both in a Monte Carlo simulation and variationally using 
either wavelets and completely local fluctuations. 



transition, indicating the presence of fluctuations at all scales. In contrary, fluctua- 
tions are strong but without scaling in the disordered phase where only short-range 
correlations persist. In the low-temperature phase, the scaling is that of the massive 
Gaussian model. 



6 Conclusions 

We have shown how wavelets can be used to describe fluctuations in a statistical 
field theory. Using a rather small number of variational parameters, an approximate 
description of fluctuations in the Gaussian, the XY, and the Landau- Ginzburg model 
has been obtained. The wavelet coefficients compactly express the scaling features of 
the field theory in only a few numbers. Comparison to Monte Carlo results leads us 
to expect that wavelets are a sufficiently efficient basis to merit further investigation 
of their use in renormalization group transformations. 

As our results are only semiquantitative, we have refrained from calculating criti- 
cal exponents directly. Since the results are arrived at numerically, it is very difficult 
to extract critical exponents by differentiating thermodynamic quantities. Once an 
implementation of the renormalization group in wavelet space is in place, they can 
be extracted much easier. 

Such renormalization group transformations can be implemented by considering 
the wavelet coefficients as the basic dynamical variables of the system and then 
integrating them out scale by scale. The wavelet transform has the advantage of 
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Figure 6: Wavelet coefficient fluctuations in the Landau- Ginzburg model, calculated 
variationally and in a Monte Carlo simulation. 



having a hierarchy built-in, unlike lattices where the hierarchy is imposed ad hoc. The 
renormalization group transformation can be made either in (wavelet) perturbation 
theory, using a variational principle, or directly when restricting the model to a 
hierarchical model []12| . Work in these directions is in progress. 

We would like to thank T. Biro, P. Carruthers, M. Greiner, and G. Sawidy for 
fruitful discussions during the progress of this project, and the German National 
Scholarship Foundation for its support. 
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